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Remarks on backward uniqueness of parabolic equations and 
incompressible Navier-Stokes well-posedness 

By F. Lam 

We explain why the theory of Escauriaza, Seregin and Sverak (Russian Math. 
Surveys, 2003) on potential finite time singularity in Navier-Stokes solutions must 
be largely misapprehended. It is found that the proofs of the backward uniqueness 
theorem for parabolic equations contain technical errors. The stated validity of 
a theorem for vorticity is established on ill-informed analyses as the solenoidal 
constraint is not taken into account. There are many cases where parabolic scalings 
are erroneously applied. We briefly discuss a number of related issues. 
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1. Background 

In the Eulerian description of the motion of an incompressible, homogeneous New¬ 
tonian fluid, the momentum and the continuity equations are 

i9u/clt -I- (u.V)u = V.u = 0, (1.1) 

where the velocity vector u = u(x,<) has the components (ui,U 2 ,U 3 ), the scalar 
quantity p = p(x, t) is the pressure, the space variable is denoted by x = (xi, X2, X3), 
and A is the Laplacian. The kinematic viscosity is v — fi/p, where p and p are the 
density and viscosity respectively. These equations are known as the Navier-Stokes 
equations. 

We seek the solution of dm as an initial-boundary value problem in the space- 
time domain denoted by HoX [0, t<T] where 0<T<oo, and Hg, in the present study, 
denotes the whole space or the half-space . We are interested in flow evolution 
from given data of finite energy. The initial condition is given by 

u(x,t=0) = Uo(x) e C^{Q.o), (1.2) 

and V.Uq = 0, in the case of the half-space, the no-slip boundary condition 

u(x, t)=0 Vz = 0. (1.3) 

Taking divergence of the momentum equation in dm and making use of the 
continuity, we obtain a Poisson equation for the pressure 

a.) 

l,j=l 

We would like to emphasise the fact that, at every instant of time, equation dH 
does not describe the dynamic evolution of the pressure as it contains no time 
information. 
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A physical (vector) quantity known as vorticity, a; = V x u, is related to the 
dynamics of angular momentum in flow motions. In view of the vector identity 
Vx(VA) = 0, the concept of vorticity ables us to consider flow motions with three 
unknowns instead of four. Vorticity evolution is governed by the vorticity equation 

— i/Aw = (w.V)u — (u.V)a;. (1.5) 

The vorticity field inherits the velocity solenoidal property of dm 

V.w = 0. (1.6) 

This is a direct consequence of the vector identity V.(VxA) = 0. For the Cauchy 
problem, the initial vorticity is specified as 

a;(x, 0) = Vxuo(x) = a;o(x) S (^“(ilo)- (1-7) 

In the case of the half-space, the value of vorticity on z = 0 is unknown a priori 
and must be determined as part of the solution. 

The velocity can be recovered from the Biot-Savart relation 

It is instructive to notice that the relation as well as the pressure given by dH 
illustrate the mutual global dependence of the flow variables at every instant of 
time. Under the hypothesis of incompressibility, any space-time perturbation in the 
vorticity will instantaneously propagate to every part of and thus modify the 
velocity (and the pressure) which in turn may generate vorticity gradients. Hence 
a local energy, the energy of fluid over a finite space, can only exist momentarily in 
flows of bounded energy except in vorticity-free motions. Nevertheless, the dynamic 
properties of the velocity and the pressure can be derived from the momentum 
equation in (HH) or from its time derivatives. 

By observation, governing equations dm and dm are invariant under time 
translation t ^ t + ts for finite constant ts which, clearly, refers to starting in¬ 
stant. The time symmetry is just Noether’s invariance theorem for energy conser¬ 
vation. Because the equations describe irreversible process of diffusion in accor¬ 
dance with the second law of thermodynamics, fluid dynamics in backward time 
[— To, 0], 0<To<oo must be investigated in terms of time-translated equations of 
motion over interval [0,To]. Fundamentally, evaluations of time-shifted solutions, 
either weak or ancient, do not shed new light on regularity. 

Given a positive parameter A, one can perform algebraic manipulations of dm 
and dm by means of transformations x —)> (x/A) and t —?> (t/A^). The dimensional 
analysis shows that if a triplet (u,p, w) solves the Navier-Stokes system, so does 
{ux,p\,uj\), where 

ua(x, t) = A“^u (A“^x, X~^t), 

px(x,t) = X~^p{X~^x,X~^t), (1.9) 

LVxi^jt) = X~‘^uj (A“^x, A“^t), 

as fluid properties p and p, are assumed homogeneous. Similar scalings apply to 
the heat equation dtQ — kAQ = 0, where Q may be interpreted as temperature. 
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Thus Qa(x, t) = Q{X~^x, The parabolic scalings are relevant in transforming 

one set of continuous triplet solution into a dilated one. The space-time dilation 
or contraction is limited by the fundamental dimensions of matters; the scaling 
parameter A cannot be arbitrarily assigned. Formally, triplet (ua,pa,wa) —>■ 0 as 
A —>■ 00 for bounded solution (u,p, w). The limit is consistent with the continuum 
hypothesis. If a triplet (u,p,a;) is shown to exist only in a space of equivalence 
classes of functions, for instance L'^{Q.q), 1 < q < oo, it makes no sense to apply the 
scalings to this triplet at any specific space-time location. 

Attempts have been made on the case of invariant norm ||u|| 2 , 3 (r 3 )/,oo[q 7 ^^). 
Specifically, Escauriaza et al (20036) claimed that the local-in-time regularity un¬ 
der the invariant scalings could be settled by means of a backward uniqueness of 
parabolic equations in conjunction with certain Carleman inequalities. The regular¬ 
ity interval is denoted by [0,Ts), where 0 < Eg < 00 . Should a singularity develop 
at time Tg from initial data (Ea, the solution would blow up at all of possible flow 
scales or 

( 1 . 10 ) 


lim sup||u(-,t)|| 7 , 3 (R 3 ) 

t ^ 1 s 


00 . 


Our priority is to draw our attention to the analytical methods involved as any 
new techniques of attacking the regularity problem ought to attract wider interest. 


2. Backward uniqueness for parabolic equations 

Escauriaza et al (2003a) intended to show a backward uniqueness theorem for the 
solution of = 0 , where the parabolic operator 


= i9t — A -I- b{x, t)X/x + c{x, t). 


( 2 . 1 ) 


It is assumed that the solution u is bounded and decays at infinity. The functions 
6 and c are measurable and bounded satisfying certain regularity and boundedness 
constraints. Denote Q.b = , where the hyper-sphere Bn = {x G M" : \x\ < 

R}. Suppose that u = u{x,t) is a bounded solution of = 0 in the parabolic 
domain Q.t = kls x [0,E > 0]. If u{-,T) = 0 in Db, then u = 0 in Q.t. Of course 
we are mainly interested in the case of n = 3. 

The proof of the backward uniqueness depends crucially on an inequality of 
Carleman type, i.e. equation (1.5) of Escauriaza et al (2003a). Eor definiteness, we 
write it down in full: There is a constant ao = ao{R,n) such that the inequalities 


\\e^u\\L2(aT) + 


< \\e^{dtu + Au)\\L2(nT) + Vm(-, T)\\L2(^n) 


( 2 . 2 ) 


hold for all a > ao and u G C^{Q,t) satisfying u{-,0) = 0 , and in the exponential 
weight function Z = a{T — f)(|a:| — R) + \x\'^. 

On pages 153-155 for the proof of Lemma 4, the last displayed equation reads: 


^ exp (^ar — -I- exp ^ — aaeR^ , 


where the domain Cl' = Br\B(i^iQa)R x [0,e > 0]. The time interval is sufficiently 
small. It was claimed that u = 0 in £iBx[0,e] in the limit of first r —>■ -too, then 
a —>■ -too and finally a —>• O’*'. 
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In fact, the limit is ambiguous because 

lim exp (r(a — r)') —>• exp f(+oo) x (— 00 )') or lim ^ 22. 

r ^+00 \ ^ ' / \ J r->+cx3 exp(r^) 00 

for finite positive a. If we let both r —>■ +00 and a —>■ +00 as these two parameters 
are independent of each other, we obtain limr,a^+oo —t exp(oo — 00 ). All of these 
limits are inderminant in the analysis of the extended reals. 

Furthermore, suppose that r remains finite while the parameter a is allowed to 
increase arbitrarily to infinity. Evidently, the second exponential function on the 
right vanishes for finite positive a and e while the first exponential tends to +c». 

There are similar technical inconsistencies in the decay limit on the last term 
on the right-hand side of equation (2.11). The term is bounded by 

a~^ exp(-r2) ||u||l=o(S4,\b2,,)x[o.i] ^ > 0. 

This bound tends to zero as the exponential function decays in the limit of r —> c», 
only if a is non-zero for fc > 1 so as to avoid indeterminancy c» x 0. However, the 
space-time parabolic cylinder 

lim (i?4r\i3r) X [0,1] —> 0. 

r—^oc 

Because the limit for ||m||l 2 does not exist as r^oo regardless whether param¬ 
eter a tends to 00 or remains finite, solution u in Lemma 4 does not necessarily 
vanish in the parabolic cylinder. We assert that the proof of the backward unique¬ 
ness for parabolic equations (Theorem 1) in the punched cylinder is incomplete and 
unjustified. 

3. Weight function in anisotropic Carleman inequality 

To apply the backward uniqueness to the Navier-Stokes equations, a Carleman 
inequality, with a revised weight function compared to (12.2p above, is introduced in 
“the half-space” R". Following Escauriaza et al (20036), we state their proposition 
6.2 or equation (6.12) as follows: Let 

4>{x,t) = 4>^^\x,t) -\- 4>^‘^\x,t), 

where 

\F\‘^ ™2c( 

(t>^^\x,t) = - — ,((>(2)(x,t) =a(l-t)-^, 

x' = {xi,X 2 ,--- ,Xn-i) SO that X = (x',Xn), and Bn = (0,0, ••• ,0,1). Then, for 
any compactly-supported function 

liG C'“(£li;R’"), 

where Hi = (M" -|-e„) x (0,1), and for any number a > ao(a), the following inequal¬ 
ity is valid: 

/ exp (2(()) (a|Mp/t^-|-|Vu|/t)da;dt 
Jai 

< c* / t^ exp (2 (/))\dtu -\- Awpdxdt. 

Jqi 
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Here c* = c*(q;) is a positive constant and a G (1/2,1) is fixed. 

The symbol a here should not be confused with the one as discussed in the pre¬ 
vious section. We emphasise the fact that the specific nature of the weight function 
is essential; the entire function exp(2(/)) in the above inequality must be bounded 
in order for the Carleman inequality to be meaningful. In their proof of the propo¬ 
sition, function A 2 was defined by 

A 2 = _ |v,/,(2)|2/^^ 


and, it was shown that 

A 2 > a -—xl--^(2a - 1) (a - 2a(2a - 2)(2a - 3)) ^ > 0. 

Since A 2 depends on the weight and hence on the number a (in fact it is a 
free positive parameter), then the following inequalities must hold 

( d \ ^ 

-- ) > K >0 

for all o > 2. As o —>■ -boo, the left-hand reduces to A 2 -b 1°° which is an inde¬ 
terminant form. Consequently, integral I in Proposition (6.1) does not exist in the 
limit. The contradiction suggests that the (arbitrary) parameter a in the Carleman 
inequality must he finite so that, at least, the weight needs to be bounded in any 
non empty open set C M" while time t is less than unity. 

One of the implications of a finiteness is that the proof of Lemma 5.3 is techni¬ 
cally flawed. The last displayed equation in the proof is 

J exp(2a())B)(|M|^-b|Vu|^)(Ay, A^(s-l/ 2 ))(s 2 /„)^exp 

/ / Q \ 2q: \ 

<<iexp(-2“-‘(- + 2)__ ”) 

where d = d{X, c*) is a bounded function, and A G (0, l/-\/T2], integration domain S 
denotes the parabolic cylinder (R”-b (3/A-bl)e„) x (1/2, 1), and the cut-off function 


X = x{yn,s) 


1 if iVn, s)eW, 
0 if (2 /„,s)^IT, 


and 


W = 



yn >1, 1/2 < s < 1, a(j)B{yn,s) < 




It was asserted that u{x,t) = 0 for all x G M" and for all t G (0,1/24) because, 
in the limit a —>■ -boo or effectively -boo, the exponential function rendered 

the right-hand side to zero. We believe that this operation of taking the limit is 
misleading; we have just shown that the Carleman inequality (6.12) has no meaning 
when (/ —>■ 00 . On the other hand, it is not clear how the analyses in the proofs 
of Lemmas 5.2-5.4 may determine the required vanishing condition u{x,t) = 0 for 
any bounded a. Certainly, the backward uniqueness for the heat equation in the 
half-space as stated in Theorem 5.1 is unattainable for unbounded weight fi. 
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4. Incompressibility and parabolic scalings 

In the problem of Navier-Stokes regularity in three space dimensions, it is vorticity 
equation CH) which may be considered as a heat equation as Sluj = 0. As the vortic¬ 
ity field is solenoidal (cf. (11.611 b the parabolic solution in any backward uniqueness 
and any Carleman inequality (the function denoted by u) must be divergence-free. 

In the proof of proposition 6.1 of Escauriaza et al 20036, the scaled function v 
is defined by v{x,t) = (exp((())u)(a;, t). The weight depends on space and time: 

^(x, t) = -|xp/(8t) - (a -f 1) log(6.(t)), 

where a > 0 is a real number, and the time-dependent function h = t exp((l — t) /3) 
for t G (0,2). Clearly .(j) ^ 0 unless x = 0 and thus neither u nor v can be 
solenoidal. (We have used the tilde symbol to distinguish the weight from the one 
defined in proposition 6.2.) By the same token, the weight function of proposition 
6.2 as written above V.cj){x) ^ 0 except at the space origin for every instant t. As a 
result, the scaled function u(x, t) = {exp{(j))u){x, t) used throughout equations (6.2)- 
(6.6) is not solenoidal because a solution u of the parabolic equation = 0 is not 
necessarily divergence-free under bound (5.1) and decay (5.3). The only exception 
is the trivial solution ui = u = v = 0. Our discussion implies that the conclusions 
drawn from the formal application of the backward uniqueness Theorem 5.1 to the 
Navier-Stokes vorticity, as stated in their §5, are simply false. 

It was stated that the pressure was bounded (cf. (3.4)) 

p € L°°{0,T). 

This is incorrect. The a priori bound obtainable by time integration of the pressure 
solution from (3.3) is no more than the summability: 

pG l3/2(K3) L^(0,T). 

Evidently, the bound of Caffarelli-Kohn-Nirenberg, p G L®/‘^(0,T), is an 

improvement. In view of embedding (time: 1 < 5/4 < 3/2), the pressure does not 
possess all of the regularity as required by Theorem l.f. 

There are a number of occasions where the meaning of the space-time norms 
(denoted by LsjiQr)) causes confusion. To be specific, consider the velocity in §3 
(cf. (3.2)): V G L^(]R^) L'^{0,T). In order to show Holder continuity, v was scaled: 
v{x,t) = Rv{x,i), where x = xq -I- Rx, and t = to + (i?>0). Evidently, the 

velocity was treated as functions. So let us consider 

v{x,t) ^ t £ > 0. 

Clearly the L^-norm of v exists over i3ix(0,1) but v is unbounded at x = —R/xq. 
Recall the fact that regularity of the weak solutions has not yet been proven at this 
stage of the analysis. In general, if u G L'^(i?ix(0,1]), 1 < g < oo, the singularity 
of counter example, v{x,t) ^ (x 6)^^“^^/^ at the space origin, demonstrates the 
fact that exercising parabolic scalings im to the weak solutions is unjustifiable. 
Regrettably, the proof of Theorem 1.3 is misguided. 

In theory of parabolic equations, fundamental solutions for operator m are 
established only for time t > 0 due to irreversibility (for adjoint in reversed time). 
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In homogeneous media filling the fundamental solution does not suggest an 
emergency-stop mechanism for viscous diffusion, which is in essence energy-dissipative 
on fine-structures in unrestrained space over time-scales substantially longer than 
macroscopic dynamics. The specification of initial data (HH), such as size and de¬ 
cay, is paramount in initial value problem (HI). Guided by physics and experiment, 
it is reasonable to envisage a reality where a proposition of backward uniqueness 
independent of problem’s initial (and boundary) conditions would not exist in all 
likelihood. Basically, dilating Navier-Stokes solutions into negative time zone by 
parabolic scalings im violates the second law of thermodynamics. (A scaled solu¬ 
tion that disrespects the law may be bounded.) 

5. Outlook 

We find that the proofs for the local regularity and potential blow-up stated by 
Escauriaza, Seregin and Sverak (2003a; 20036) contain analytical mistakes as well 
as misconceptions of serious nature. Claims and applications made on the basis of 
the proposed problem-solving scheme require critical reappraisals. 

Any incompressible fluid motion involves velocity, pressure as well as vorticity; 
the aggregate flow field is scale-sensitive in norm L^. Equally, we may as well 
consider the vorticity equation in Lx • There does not exist one single definitive 
scale in which triplet (u,p, w) is invariant during flow evolution. On the other hand, 
the energy identity of regular Navier-Stokes solutions, 

1 11 11 2 II 11 2 ^11 11 ^ 

2 II'^(^)|Il2(r 3) ||^u(f)||^2^jj3jdf = 2 II'^o|Il2(r3)’ 

is truly scale-independent. Both the velocity and its Jacobian are taken into account, 
and, implicitly, the pressure or the continuity is also effected. The conservation law 
does not demand any scale-dependency on initial data dni). By virtue of the Biot- 
Savart relation m or the pressure gradient from dOD, triplet (u,p, w) is essentially 
non-local. The solutions are known globally regular; no singularity is possible on 
any local scales. Experience shows that turbulent motions decay in time. The last 
stage of the decaying flow may be modelled approximately by diffusion equation, 
i9tu— i/Au = 0 that has only constant-scale solutions. Thus parabolic scalings (jl.9ll 
do not play a significant role in understanding turbulence whose spatio-temporal 
structures are theoretically tractable by a proper treatment of the non-linearity in 
the governing equations. The finite-time singularity scenario il.KM) is illusive. 


17 September 2015 
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